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ABSTRACT In this paper, a control chart for variance using the repetitive sampling under the neutrosophic 
statistics is proposed. The proposed control chart can be applied for monitoring the variability in the 
process when some observations/parameters are unclear and indeterminate. The operational procedure of 
the proposed control chart and the necessary measures are given for the practical use of the proposed control 
chart. The control chart parameters are determined through the algorithm under neutrosophic statistical 
interval method. The efficiency of the proposed control chart is compared with the existing control chart 
in neutrosophic average run length. From the simulation study and a real example, it is concluded that the 
proposed control chart is more efficient in detecting a shift in the process and more adequate, effective, 


informative, and flexible under the uncertainty environment. 


INDEX TERMS Neutrosophic statistical interval method, neutrosophic average run length, variance, 


repetitive sampling, uncertainty, control chart. 


I. INTRODUCTION 

To meet the high standard of quality of the product, it is 
necessary to monitor the variation in the process. The pres- 
ence of the variation in the process may shift the process 
from the set target. According to [1] “High quality in a 
product can be only achieved if the manufacturing process 
meets the given specifications limits. To the manufacturer to 
produce the defect-free product, the variation in the process 
should be controlled”. A control chart is an important tool, 
which has been commonly used for the monitoring variation 
in the process. The control chart is a graph having lower 
control limit (LCL), central limit (CL) and the upper control 
limit (UCL). The process is said to be out-of-control if the 
plotting statistic is beyond LCL or UCL. The late indication 
about the shift in the process causes the production of non- 
conforming items. The control chart immediately indicates 
when the process is shifted from the specified target. There- 
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fore, the control charts are designed to help the operators to 
fix the issue in the manufacturing process before the deeper 
problem in the process. The Shewhart S 2 control chart has 
been commonly applied for the monitoring of the variation in 
the process. The S* chart using the double sampling was pro- 
posed by [2]. Reference [3] proposed the ARL-unbiased con- 
trol chart. Reference [4] designed the efficient control chart to 
monitor the variation in the process. Reference [5] worked on 
the variance control chart using combined sampling schemes. 
Reference [6] compared the efficiency of interquartile range 
chart with existing control charts. For more details, the reader 
may refer to [7]-[9]. 

The control charts are usually designed using the single 
sampling scheme. In this type control charts, a decision about 
the state of the control chart is decided based on sample 
information from the single sample selected from the pro- 
cess. The repetitive sampling (RS) scheme was originally 
proposed by [10] which is applied when the experimenter is 
indecisive about the state of the process at the first sample. 
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The repetitive sampling based control charts are found to be 
more efficient in average run length (ARL) than the control 
chart using the single sampling. Reference [11] designed 
control chart for process capability using the RS. Reference 
[12] worked on t-chart using this sampling scheme. Refer- 
ence [13] deigned sign chart using the RS. Reference [14] 
proposed a hybrid type control chart using the RS. Reference 
[15] proposed the S? control chart using the RS. More details 
on the RS and ranked sampling control charts can be read 
in [16]-[21]. 

The fuzzy control charts, which are more sensitive than 
the Shewhart control, are applied when parameters are 
uncertain [22]. Reference [23] “observations include human 
judgments, and evaluations and decisions, a continuous 
random variable of a production process should include 
the variability caused by human subjectivity or measure- 
ment devices, or environmental conditions. These variabil- 
ity causes create vagueness in the measurement system”. 
Reference [23] proposed the dispersion control chart using 
the fuzzy logic. Reference [23] proposed a variance control 
chart using a-cut. [24] designed control chart under uncer- 
tainty. The applications of the control chart using the fuzzy 
approach can be seen in [25]. More details can be seen 
in [26]-[29]. 

Reference [30] mentioned that the neutrosophic logic 
which deals with indeterminacy is the generalization of the 
fuzzy logic. Reference [31] used the neutrosophic logic to 
develop the idea of the neutrosophic statistics (NS). The NS 
can be applied for the analysis of the data, which is drawn 
from pupation or sample having uncertain observations or 
the parameters. The NS is the extension of the classical 
statistics and found to be more effective to analyze the 
data under the uncertainty. The NS provides the analysis 
results in indeterminacy interval rather than the determined 
values. Therefore, the NS is more adequate and effective 
to a method to apply under the uncertainty environment. 
References [32] and [33] suggested that NS is a more effec- 
tive method than the classical statistics in rock measuring 
issues. References [34] and [35] originally introduced the 
NS in the area of acceptance sampling plans. Reference [36] 
introduced the NS in the area of control chart first time. 
Reference [1] proposed the variance control chart using 
the NS. 

By exploring the literature and best of our knowledge, 
there is no work on the S? control chart using the repetitive 
sampling under the NS. In this paper, we will design a vari- 
ance control chart using the neutrosophic statistical interval 
method (NSIM). We expect that the proposed control chart 
will be more efficient in detecting a shift in the process 
and more adequate, effective, informative and flexible in the 
uncertainty environment. 


Il. PREMIERES 


Let Xni € {XL, Xu}, i = 1,2,3,..., my be a neutrosophic 
random variable (NRV) consisting of lower value Xz and 
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FIGURE 1. The proposed control chart for the simulation data. 


upper value Xy. Suppose also that uy = pn Xyi/ NN; 
uy € {uL nu} and og = DN, Æi — uN)’ / (Ny — 1); 
of € {o?, aa} represent the neutrosophic population mean 
and variance corresponding to Xz and Xy, respectively. The 
neutrosophic estimators of the population mean is Xy = 
paar Xni/tn; XN € XL, Xu} and population variance 
; =. V2 

is Sy = Di (Xm — Xv)" / ww — Ds SH E {S SG}. 
From [1] and [31], Sy E€ Is, S2} has approximate neu- 
trosophic normal distribution with neutrosophic mean of € 
{o?, of} and neutrosophic variance of) / (ny — 1). Also, 
(ny — 1) Si /oqs Sh € {S7,S?},09 € {o}, of} is mod- 
elled by the neutrosophic Chi-Square distribution (x aa Xo € 
Fen xa) having a neutrosophic degree of freedom (ndf) 
vy = ny — l; ny € {nL, ny}. Also, Gy € {GL, Gy} presents 
the neutrosophic cumulative distribution function (ncdf), see 
for example, [1]. 


Ill. DESIGN OF THE PROPOSED CHART 

Based on the given information, we propose the follow- 
ing control chart for the repetitive sampling under the 
NSIM. 

Step-1: Compute S? based on NRV of size ny. 

Step-2: Declare the process is an in-control state if 
LCLy2 < Se < UCLy2; where LCLy2 € {LCLz2, LCLy2} 
and UCLy2 € {UCL 12, UCLy2} are neutrosophic upper 
control limits. 
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FIGURE 2. The existing chart for the simulated data. 


Step-3: The process is said to be out-of-control if S3 > 
UCLyı or SẸ < UCLy1; where LCLy1 € {LCL11, LCLy1} 
and UCLyı € {UCLz1, UCLyı} are neutrosophic lower 
control limits, otherwise repeat the process. 

The proposed control consists of four control limits with 
two-neutrosophic control limits coefficient under the NSIM. 
The proposed chart shows the control limits in indeterminacy 
interval. Under this complex system, the operational proce- 
dure of the proposed control chart is the same as in double 
sampling and repetitive sampling under classical statistics. 
Also, the proposed control chart is the generalization of the 
control chart proposed by [1] under the NSIM and [15] 
when no certainty in the observations or in the parameters. 
The neutrosophic control limits LCLN; € {LCLz1, LCLy1}, 
UCLyı € {UCL11, UCLy1}, LCLn2 € {LCL:12, LCLy2}, 
UCLy2 € {UCLr2, UCLy2} of the proposed control chart 
are given as 





LCLyı = on = wa (02) / (nw = 1); 


oh efoz.og}, kvie tknku) ®© 





UCLyy = of + kyiy 2 (02) / (nw — 1); 


ok e {oz. 03}. kyi € {kr1, kui} (2) 





LCLw2 = 02 — kvay/2 (o2)° / (aw — 1); 
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Gee {o?. og} » kn2 € {kr2, ku2} (3) 





2 
UCLy2 = of + kwoy/2 (of) / 0n — 1); 
oh € [o2 og} » kyo € {kr2, ku2} (4) 


Note here that ky; € {kz1, ky1} and ky2 € {kz2, ky2} are the 
control limits coefficients. The probability when the process 
is out-of-control using the single sampling under the NSIM 
is given by 


P(S} > UCL) 
—1)UCL 
Sisi (= ) “) 





2 
ON 


| 2: 
= 1—Gy (1 <9 (140 a): 


oh € [o2 of} , kyi € {kr1, kui}, ny € {nz, nu} (5) 
or 


P (Sà < LCLy1) 


(= = ucim) 
= N ~n 


ON 
= Gy ((ny - D (1 = kvıV2/ 0N = D) ); 
oe (oz, of} » kyi € {kz1, kui}, ny € {nL, nu} (6) 


Suppose that the process is declared to be an in-control state at 
of € {o? ; og}. Therefore, the probability of out-of-control 
for the single sampling under the NSIM can be written as 


follows 


0 = = a 
Pauni = 1—Gn (ox 1) (1440 Ge 5)) 
+G 1 l-k 2 ; 
N { (nN ) N1 Gy aD : 


on e {of o8}, kni € {kii, kui}, nn € {nL, nu} 
(7) 


The neutrosophic probability of indecision at the first sample 
is given by 





pov = P (UCLy2 < S2 < UCLNi lox) 


+P(LCLy1 < Sy < LCLy2lo%) 


one fop ob}, kvietku) ® 


The Eq. (8) can be written as follows 


0 _ E | 2 
P reptN = Gn (o 1) ( + kn Gn- D — 5)) 
G D{1+k se 
— Gy | (nn — 1) | 1 +k m-i) 
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TABLE 1. The values of NARL when ny € {3, 4}. 




































































TABLE 2. The values of NARL when ny € {4, 6}. 
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ron=300 ron=370 
Ny [3,4] [3,4] 
ky [(5.8737, 3.9607),( 5.2171, [(6.0840, 5.6889),( 5.4234, 
5.1036)] 4.5386)] 
C 
1 [300.93, 300.15] [371.99, 374.31] 
1.1 [166.72, 159.31] [203.55, 193.75] 
1.2 [101.98, 94.51] [123.48, 112.45] 
1.3 [67.31, 61.04] [81.06, 71.23] 
14 [47.16, 42.13] [56.61, 48.31] 
15 [34.65, 30.66] [41.53,34.59] 
1.6 [26.47, 23.28] [31.71, 25.89] 
1.7 [20.88, 18.30] [25.02, 20.08] 
18 (16.92, 14.81] (20.29, 16.05] 
1.9 (14.02, 12.28] (16.83, 13.16] 
2 [11.85, 10.39] (14.24, 11.02] 
3 [4.18, 3.78] [5.05, 3.76] 
4 [2.58, 2.39] [3.08, 2.32] 
ron=300 ron=370 
Ny [4,6] [4,6] 
ky [(5.2314, 3.5838),( [(5.4133,4.5360),(4.7725, 
4.6238,4.5484)] 4.6760)] 
Cc 
1 [300.17, 309.25] [370.21,374.86] 
1.1 [157.81, 150.29] [191.84,178.40] 
1.2 [92.64, 83.35] [111.45,97.23] 
1.3 59.16, 51.10 [70.65,58.74] 
1.4 [40.36, 33.87] [47.95,38.45] 
1.5 [29.02, 23.88] [34.36,26.81] 
1.6 21.78, 17.69 [25.72,19.68] 
1.7 16.94, 13.64 [19.96,15.05] 
1.8 13.56, 10.88] [15.97,11.91] 
1.9 [11.13, 8.92] [13.09,9.70] 
2 [9.33, 7.48] [10.97,8.09] 
3 [3.24, 2.71] [3.75,2.83] 
4 [2.04, 1.78] [2.32,1.83] 














1 


k 2 
NIV (ay — 1) 





a T 


k : : 
N2 aD)’ 


2 2 ~2 





€ {kr1, kui}, 


ny € {nL, ny} (9) 


Finally, the probability of out-of-control for the repetitive 
sampling under the NSIM is given by 


VOLUME 7, 2019 


M. Aslam: Control Chart for Variance Using RS Under Neutrosophic Statistical Interval System 


IEEE Access 





TABLE 3. The values of NARL when ny € {9, 10}. 



























































ron=300 ron=370 
Ny [9,10] [9,10] 
ky [(4.1807, 3.4441),( 4.0979, | [(4.3148, 4.2484),( 4.2244, 
3.4762)] 3.3250)] 
C 
1 [301.71, 305.75] [370.62,370.46] 
1.1 [131.97, 129.72] [159.14,153.15] 
13 [67.48, 64.88] [80.35,74.88] 
1.3 [38.82, 36.73] [45.85,41.54] 
1.4 [24.46, 22.89] [28.76,25.41] 
1.5 [16.57, 15.38] [19.43,16.79] 
1.6 [11.89, 10.97] [13.94,11.80] 
1.7 [8.94, 8.23] [10.48,8.73] 
1.8 [6.99, 6.42] [8.20,6.74] 
1.9 [5.65, 5.19] [6.63,5.38] 
2 [4.69, 4.31] [5.51,4.43] 
[1.76, 1.65] [2.01,1.64] 
4 [1.28, 1.23] [1.40,1.22] 











p? 

0 _ outN| . 0 0 0 

Poun = 1 po , Poun € eae Pou | (10) 
~*~ reptN 


The neutrosophic average run length (NARL) pro- 
vides the indeterminacy interval in which first out-of- 

: 2 2 421; 
control is expected. The NARL at oy € {ar og} is ispiven by 
given by 


NARLon = 





g; ARLon € {ARLoL, ARLou} (11) 


We suppose that due to some external variations, the process 


: : 2 = Da 329 2 2 

is shifted at new value On, = CON; On, €E ae oi}. 
where c shows the shift constant. Therefore, the probabil- 
ity of out-of-control at On = cof; Ge, E€ lof oki} 





(ny — 1) 2 : 
+ Gy z ( T] 


ki € {kr1, kui}, 


ny € {nL, nu} (13) 


The neutrosophic probability of indecision at the first sample 





1 Co (ny — 1) 2 
Po Preptn = 2 (o ( A ( +kyı N) 5)) 
(ny — 1) | 2 ; 
— Gy (=> ( + ky ar- 5))) ; 


kn2 € {kr2, ku2}, 


ny € {nL, ny} (14) 


for the single sampling under the NSIM can be written as 
follows 


Phun = P (Sh = UCLyilog) +P (Si < LCLwilog) ; 
oN € {oz of: (12) 


or 


ee (ny — 1) 2 
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Finally, the probability of out-of-control for the repetitive 
sampling under the NSIM is given by 


p! 

1 = outN1 . 1 1 1 

outN 7T 1 pP! , Poun € | Pic Prat} (15) 
~~ reptN 


The NARL at oži € lof oka} is given by 





NARLin = ARLin € {ARLiL, ARLiv} (16) 


1 ; 
outN 
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TABLE 4. The comparison in NARL when ny € {4, 6} and ARLoy € {300, 300}. 
























































c Proposed Chart Existing Chart 
NARL NARL 

1 [300.17, 309.25] [300.03,490.72] 
1.1 [157.81, 150.29] [167.75,236.62] 
1.2 [92.64, 83.35] [103.63,130.01] 
1.3 [59.16, 51.10] [69.11,78.91] 
1.4 [40.36, 33.87] [48.93,51.76] 
1.5 [29.02, 23.88] [36.34,36.10] 
1.6 [21.78, 17.69] [28.05,26.46] 
1.7 [16.94, 13.64] [22.35,20.20] 
1.8 [13.56, 10.88] [18.28,15.94] 
1.9 [11.13, 8.92] [15.29,12.93] 
2 [9.33, 7.48] [13.03,10.74] 
3 [3.24, 2.71] [4.84,3.58] 

4 [2.049, 1.78] [3.02,2.21] 








Suppose that NARLon be specified by roy. The values 
of NARL\n for various combinations of ny € {nz, ny}, 
ARLon € {ARLo, ARLou} and c are presented in 
Tables 1-3. From Tables 1-3, we note the decreasing trend in 
the indeterminacy interval of NARL;yn when c changes from 
1 to 4. The indeterminacy interval of NARL jn is increasing 
when ARLon changes from 300 to 370. The indeterminacy 
interval of NARL n is decreasing when ny € {nz,ny} 
increases. 

The following algorithm is used under the NSIM to deter- 
mine kyı € {kzi,kyi},knw2 € {kr2, ku2} and ARLin € 
{ARL;L, ARLiv}. 

Step-1: Pre-fix the indeterminacy interval of ny € 
{nL, ny} and shift value c. 

Step-2: Specify the suitable ranges for ky; € {kz1, ku1}, 
ky2 € {kr2, ky2}, where kr < ky and k2 < ky2. Deter- 
mine the indeterminacy interval of ky, € {kz1, kui}, ky2 € 
{kr2, ky2} such that NARLon > ron. 

Step-3: During, the simulation, we note that several combi- 
nations of ky, € {kr1, ky1},kw2 € {kr2, ky2} are available 
that satisfy the condition NARLon=ron. 

Step-4: choose that values of kyı € {kz1,kui},kw2 € 
{kr2, ku2} where NARLon is same or very close 
to ron. 
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Step-5: Determine ARLin € {ARLi_, ARLjv} for various 
values of c. 


IV. ADVANTAGES OF THE PROPOSED CHART 

In this section, we will discuss the advantages of the proposed 
control chart over the control charts proposed by [1] under 
the NSIM. We will compare three-control chart in terms of 
NARL and using the simulation data. 


A. COMPARISON IN NARL 

We will now compare the proposed chart under the NSIM 
with chart proposed by [1] under the NSIM in NARL. A con- 
trol chart having the smaller values of NARL at the same 
parameters is said to be a more efficient control chart. Let 
R, = ARLiu — ARL; be the neutrosophic range of inde- 
terminacy interval of NARL. The smaller the values of R, 
better the control chart. The values of NARL of both control 
charts when ny € {4, 6}, ARLon € {300, 300} and ARLon € 
{370, 370} are shown in Tables 4-5. We compare both control 
chart in the rage of indeterminacy interval of NARL. From 
Tables 4-5, we note that the proposed control chart under 
NSIM has smaller values of NARL as compared to [1] control 
chart. For example, when ny € {4, 6}, ARLon € {300, 300} 
and c = 1.2. The values of indeterminacy interval of NARL 


VOLUME 7, 2019 


M. Aslam: Control Chart for Variance Using RS Under Neutrosophic Statistical Interval System 


IEEE Access 








0.025 
1 
o, 





=- S 
L 
UCL, ;=0.021628 


UCL; =0-0205335 











0.020 
1 


UCL, ,=0.018377 





UCL y2=0.016771 


. 
. i 
P 
ile 
. 


. 
. 





0.015 
| 
SSS ite 


S; 
0.010 
1 





0.005 
1 











0.000 


RERE 


LChy=—0-0038 








LOL y=+0:00283 








-0.005 
1 





Sample Number 


FIGURE 3. The proposed control chart for the real data. 


and R, of proposed chart are ARLin € {92.64, 83.35} and 
Rn = 9. On the other hand, values of indeterminacy interval 
of NARL and R, of [1] chart are ARL1n € {103.63, 130.01} 
and R, = 27. From these values, we note that the proposed 
control chart will give the indication in the shift between 
83th and 92" samples while the existing control chart will 
provide the single between 103th and 130 samples. The 
By comparing the proposed control chart with [1] control 
chart, we conclude that the proposed control chart is more 
efficient than the existing control chart in detecting a shift 
in the process. Therefore, the proposed control chart is more 
effective than the existing control chart under the NSIM. 


B. COMPARISON IN NARL 

Now, we compare the efficiency of the proposed control 
chart under the NSIM over [1] using the simulation data 
generated from the neutrosophic normal distribution with 
neutrosophic mean and variance wy € {0,0} and on E€ 
{4, 6.25}. First 20 observations are generated assuming the 
process is in-control and next 20 observations are generated 
for the shifted process when c = 1.9. Let ny € {5,5} and 
ARLon € {370, 370} for this study. The values of statistic 
Sx € 1s; S? }, are computed and plotted on Figure 1. From 
Figure 1, it can be seen in that the proposed control chart 
detects shift at the 8 sample. On Figure 2, the values of 
statistic Si € fs r, s?) are plotted on the existing control 
chart. From Figure 2, it is noted that the existing control 
chart detects the first shift in the process at the 14" sample. 
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FIGURE 4. the existing control chart for the real data. 


By comparing Figure 1 with Figure 2, it is concluded that the 
proposed chart has the ability to detect an early shift in the 
process as compared to [1]. 


V. CASE STUDY 

A well-known automobile company in Jeddah, Saudi Arabia 
is interested to apply the proposed control chart for the mon- 
itoring of engine piston rings. The diameter of engine piston 
rings is the quality of interest which is an important vari- 
able, see [37]. A small increase or decrease in the diameter 
may cause the defective product. The instrument does the 
measurement of diameter, which causes the uncertainty in 
the observations. As mentioned above, measurements done 
by the devices or by human subjectivity may cause to have 
some uncertain values in the data, see [23]. The data hav- 
ing uncertain observations cannot be monitored effectively 
using [15] control chart under classical statistics. We bor- 
rowed the piston rings from [1] and reported in Table 6 for 
the easy reference. 

The values of S? E {s?, s?) are calculated and shown 
in Table 6. The neutrosophic control limits when ny € {5, 5}, 
ARLon € {370,370}, kyı € {4.8582, 4.9535}, ky € 
{3.5665, 3.9374} and of € {0.0035, 0.0034} are computed 
as follows 





LCL = o2 = kriy2 (o2)? / (ny — 1) = —0.0028 
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TABLE 5. The comparison in NARL when ny € {4, 6} and ARLoy € {370, 370}. 


















































c Proposed chart Existing Chart 
NARL NARL 
1 [370.21,374.86] [370.33, 379.28] 
1.1 [191.84,178.40] [192.50, 180.33] 
1.2 [111.45,97.23] [112.25, 98.22] 
1.3 [70.65,58.74] [71.47, 59.32] 
1.4 [47.95,38.45] [48.73, 38.83] 
1.5 [34.36,26.81] [35.09, 27.08] 
1.6 [25.72,19.68] [26.40, 19.87] 
1.7 [19.96,15.05] [20.60, 15.20] 
1.8 [15.97,11.91] [16.55, 12.04] 
1.9 [13.09,9.70] [13.64, 9.81] 
[10.97,8.09] [11.48, 8.18] 
[3.75,2.83] [4.04, 2.86] 
[2.32,1.838] [2.51, 1.85] 

















UCL = of + iy (a2) / (ny — 1) = 0.0216 





LCLr2 = o2 = ky 2 (a2)? / (ny — 1) = 0.0004 





UCL 12 = of + kay/2 (o2 / (ny — 1) = 0.0183 





LCLy = og — kay2 (02)? / (ny — 1) = —0.0038 





UCLyi = 02 + kuy2 (o2)° / (nw — 1) = 0.0205 





LCLy2 = 03 — kuoy/2.(02)" / (aw — 1) = —0.00004 





Ulur = 0? + hay 2(02) / Gn - 1) = 0.0167 


The values of statistic Sž € {$7 , S?) are plotted in Figure 3. 
From Figure 3, we can note that several points are near the 
control limits. We also note that the 2" sample is in indeter- 
minacy interval. The existing chart by [1] in Figure 4 indicate 
the process is in control. The proposed chart for this real data 
is implemented as follows 

Step-1: Compute Sx, € [0.014772, 0.017242], ...S} € 
[0.016177, 0.013088] as in Table 6 based on NRV of size 
ny € {5, 5}. 

Step-2: Declare the process is in-control state if 
[0.0004, 0.0028] < Sy < [0.0183, 0.0216]. 

By comparing Figure 3 with Figure 4, we conclude that 
the proposed control chart indicates some shifting issues in 
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the process although the process is in control state. We can 
conclude that the proposed control chart has the ability to 
detect an early shift in the process as compared to the existing 
control chart. As the proposed control chart provides results 
in indeterminacy interval rather than the determined values, 
therefore, the proposed control chart is more adequate under 
uncertainty environment. 


VI. CONCLUDING REMARKS 

In this paper, a control chart for variance using the repetitive 
sampling under the neutrosophic statistics (NS) is proposed. 
The proposed control chart is the generalization of several 
control charts. The proposed control chart using the repetitive 
sampling is more flexible and allow to repeat the process if 
the operations are in-decision at first sample information. The 
proposed control chart provides the indeterminacy interval 
for the NARL values which is more desirable under the 
uncertainty. From the comparison, we conclude that the pro- 
posed is more sensitive in detecting the shift in the process. 
The proposed control chart can be used for the monitoring 
of process in the automobile industry, aerospace industry 
and food industry. The proposed control chart for the mon- 
itoring of big data in marine science can be extended for 
future research. The proposed chart using the cost model can 
also be considered as future work. The joint monitoring of 
mean and variance under NSIM can be considered as future 
research. 
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TABLE 6. Real example data. 











































































































Sample Observation Si 
1 [74.03,74.03] | [74.002,73.991] | [74.019,74.019] | [73.992,73.992] | [74.008,74.001] | [0.014772,0.017242] 
2 [73.995,73.995] | [73.992,74.003] | [74.001,74.001] | [74.011,74.011] | [74.004,74.004] | [0.007503,0.005762] 
3 [73.988,74.017] | [74.024,74.024] | [74.021,74.021] | [74.005,74.005] | [74.002,73.995] | [0.014748,0.012116] 
4 [74.002,74.002] | [73.996,73.996] | [73.993,73.993] | [74.015,74.015] | [74.009,74.009] | [0.009083,0.009083] 
5 [73.992,73.992] | [74.007,74.007] | [74.015,74.015] | [73.989,73.989] | [74.014,73.998] | [0.012219,0.010756] 
6 [74.009,74.009] | [73.994,74.001] | [73.997,73.997] | [73.985,73.985] | [73.993,73.993] | [0.008706,0.008944] 
7 [73.995,73.998] | [74.006,74.006] | [73.994,73.994] [74,74] [74.005,74.005] | [0.005523,0.00498] 
8 [73.985,73.985] | [74.003,74.01] | [73.993,73.993] | [74.015,74.015] | [73.988,73.988] | [0.012256,0.01348] 
9 [74.008,74.005] | [73.995,73.995] | [74.009,74.009] | [74.005,74.005] | [74.004,74.004] | [0.005541,0.005177] 
10 | [73.998,73.998] [74,74] [73.99,73.99] [74.007,74.007] | [73.995,73.995] | [0.006285,0.006285] 
11 [73.994,73.998] | [73.998,73.998] | [73.994,73.994] | [73.995,73.995] | [73.99,74.001] | [0.002864,0.002775] 
12 | [74.004,74.004] [74,74.002] [74.007,74.005] [74,74.001] [73.996,73.996] | [0.004219,0.003507] 
13 | [73.983,73.993] | [74.002,74.002] | [73.998,73.998] | [73.997,73.997] | [74.012,74.005] | [0.010455,0.004637] 
14 | [74.006,74.006] | [73.967,73.985] | [73.994,73.994] [74,74] [73.984,73.996] | [0.015304,0.007759] 
15 | [74.012,74.012] | [74.014,74.012] | [73.998,73.998] | [73.999,73.999] | [74.007,74.007] | [0.007314,0.006804] 
16 [74,74] [73.984,73.984] | [74.005,74.005] | [73.998,73.998] | [73.996,73.996] | [0.007797,0.007797] 
17 | [73.994,73.994] | [74.012,74.012] | [73.986,73.986] | [74.005,74.005] | [74.007,74.007] | [0.010569,0.010569] 
18 | [74.006,74.006] | [74.01,74.011] | [74.018,74.018] | [74.003,74.003] [74,74.001] [0.006986,0.006834] 
19 | [73.984,73.984] | [74.002,74.002] | [74.003,74.003] | [74.005,74.005] | [73.997,73.997] | [0.008468,0.008468] 
20 [74,74] [74.01,74.01] | [74.013,74.009] | [74.02,74.015] | [74.003,74.003] | [0.007981,0.005941] 
21 [73.982,73.982] | [74.001,74.001] | [74.015,74.015] | [74.005,74.005] | [73.996,73.996] | [0.012153,0.012153] 
22 | [74.004,74.004] | [73.999,73.999] | [73.99,73.99] [74.006,74.006] | [74.009,74.002] | [0.007436,0.006261] 
23 [74.01,74.01] | [73.989,73.989] | [73.99,73.99] [74.009,74.005] | [74.014,74.011] | [0.011929,0.010747] 
24 | [74.015,74.011] | [74.008,74.008] | [73.993,73.993] [74,74] [74.01,74.011] | [0.008701,0.007893] 
25 | [73.982,73.982] | [73.984,73.989] | [73.995,73.995] | [74.017,74.012] | [74.013,74.01] | [0.016177,0.013088] 
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